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«) | Abstract. We present two new families of stationary solutions for equations of Bose-Fermi 

(*C) ■ mixtures with an elliptic function potential with modulus k. We also discuss particular 

cases when the quasipcriodic solutions become periodic ones. In the limit of a sinusoidal 
potential (k — ► 0) our solutions model a quasi-one dimensional quantum degenerate Bosc- 
^0 | Fermi mixture trapped in optical lattice. In the limit k —> 1 the solutions are expressed 

£H • by hyperbolic function solutions (vector solitons). Thus we are able to obtain in an unified 

way quasi-periodic and periodic waves, and solitons. The precise conditions for existence 
of every class of solutions are derived. There are indications that such waves and localized 
objects may be observed in experiments with cold quantum degenerate gases. 
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1 Introduction 



Over the last decade, the field of cold degenerate gases has been one of the most active areas 
in physics. The discovery of Bose-Einstein Condensates (BEC) in 1995 (see e.g. [U[2]) greatly 
stimulated research of ultracold dilute Boson- Fermion mixtures. This interest is driven by the 
J> ■ desire to understand strongly interacting and strongly correlated systems, with applications in 

solid-state physics, nuclear physics, astrophysics, quantum computing, and nanotechnologies. 

An important property of Bose-Fermi mixtures wherein the fermion component is dominant 
is that the mixture tends to exhibit essentially three-dimensional character even in a strongly 
elongated trap. During the last decade, great progress has been achieved in the experimental 
realization of Bose-Fermi mixtures [3j d] , in particular Bose-Fermi mixtures in one-dimensional 
lattices. Optical lattices provide a powerful tool to manipulate matter waves, in particular 
solitons. The Pauli exclusion principle results in the extension of the fermion cloud in the 
transverse direction over distances comparable to the longitudinal dimension of the excitations. 
It has been shown recently, however, that the quasi-one-dimensional situation can nevertheless 
be realized in a Bose-Fermi mixture due to strong localization of the bosonic component [5j[6]. 
With account of the effectiveness of the optical lattice in managing systems of cold atoms, their 
effect on the dynamics of Bose-Fermi mixtures is of obvious interest. Some of the aspects of 
this problem have already been explored within the framework of the mean-field approximation. 
In particular, the dynamics of the Bose-Fermi mixtures were explored from the point of view 
of designing quantum dots [8]. The localized states of Bose-Fermi mixtures with attractive 
(repulsive) Bose-Fermi interactions are viewed as a matter-wave realization of quantum dots 
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and antidots. The case of Bose-Fermi mixtures in optical lattices is investigated in detail 
and the existence of gap solitons is shown. In particular, in [5] it is obtained that the gap 
solitons can trap a number of fermionic bound-state levels inside both for repulsive and attractive 
boson-boson interactions. The time-dependent dynamical mean-field-hydrodynamic model to 
study the formation of fermionic bright solitons in a trapped degenerate Fermi gas mixed with 
a Bose-Einstein condensate in a quasi-one-dimensional cigar-shaped geometry is proposed in [9]. 
Similar model is used to study mixing-demixing in a degenerate fermion-fermion mixture in |10j . 
Modulational instability, solitons and periodic waves in a model of quantum degenerate boson- 
fermion mixtures are obtained in [11] . 

Our aim is to derive two new classes of quasi-periodic exact solutions of the time dependent 
mean field equations of Bose-Fermi mixture in one-dimensional lattice. We also study some 
limiting cases of these solutions. The paper is organized as follows. In Section [2] we give the 
basic equations. Section [3] is devoted to derivation of the first class quasi-periodic solutions 
with non-trivial phases. A system of Nf + 1 equations, which reduce quasi-periodic solutions to 
periodic are derived. In Section |4] we present second class (type B) nontrivial phase solutions. 
In Section [5] we obtain 14 classes of elliptic solutions. Section [6] is devoted to two special limits, 
to hyperbolic and trigonometric functions. In Section [7] preliminary results about the linear 
stability of solutions are given. Section [8] summarizes the main conclusions of the paper. 



2 Basic equations 

At mean field approximation we consider the following Nf + 1 coupled equations El fl~2| [TT] 

ih ~dt + 2^^ " v ^ b ~ 5bb| * 1 " gB¥p f^ b = °' (2 - x ) 

^ + 2W^-^-^'M = ^ (2-2) 

N f 

where pf = I ^{ I an d 

i=l 



2obb 2obf itib / h 

9bb = , Qbf = , a = , a s = 4 / , 

a s a s a m-p y tob^i 

obb and cjbf are the scattering lengths for s-wave collisions for boson-boson and boson-fermion 
interactions, respectively. In recent experiments |13tll4| the quantum degenerate mixtures of 40 K 
and Rb are studied where ttt-b = 87m p , me = 40m p and to± = 215 Hz. Equations (|2.ip . (|2.2p 
have been studied numerically in [7|. The formation of localized structures containing bosons 
and fermions has been reported in the particular case in which the interspecies scattering length 
(zbf is negative, which is the case of the 40 K- 87 Rb mixture. An appropriate class of periodic 
potentials to model the quasi-lD confinement produced by a standing light wave is given by [15] 



V = Vosn 2 (ax, k), 

where sn (ax, k) denotes the Jacobian elliptic sine function with elliptic modulus < k < 1. 

Experimental realization of two-component Bose-Einstein condensates have stimulated con- 
siderable attention in general [16] and in particular in the quasi- ID regime \17\ [T8] when the 
Gross-Pitaevskii equations for two interacting Bose-Einstein condensates reduce to coupled 
nonlinear Schrodinger (CNLS) equations with an external potential. In specific cases the two 
component CNLS equations can be reduced to the Manakov system |19] with an external po- 
tential. 
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Important role in analyzing these effects was played by the elliptic and periodic solutions 
of the above-mentioned equations. Such solutions for the one-component nonlinear Schrodinger 
equation are well known, see [20] and the numerous references therein. Elliptic solutions for the 
CNLS and Manakov system were derived in [2H [22], I23j . 

In the presence of external elliptic potential explicit stationary solutions for NLS were derived 
in [15 ^ 124 } [25]. These results were generalized to the n-component CNLS in [18] . For 2-component 
CNLS explicit stationary solutions are derived in |26j. 



3 Stationary solutions with non-trivial phases 

We restrict our attention to stationary solutions of these CNLS 

^ b (x, t) = q (x) exp (-i^-t + i@o(x) + Ikq) , (3.1) 



vUx,t) = qj{x) exp (-i^-t + i@j(x) + iK ,j) , (3.2) 



h 

where j = 1, . . . ,Nf, kq, kqj, are constant phases, qj and 0o, &j(x) are real-valued functions 
connected by the relation 



o (x) 



dr' f x dr' 

Jo TAx') 



<7oV) 



CojCj') j = 1, - - - , JVy being constants of integration. Substituting the ansatz (|3.ip . ()3.2[) in 
equations (|2.ip and separating the real and imaginary part we get 

Nf \ 1 

qfaoxx - 5BB<?0 - Vql - 5BF I ^2 q i Q + W °^0 = 9 S 0> ( 3 - 4 ) 



2m B u u u 2m B 



We seek solutions for q$ and q 2 , j = 1, . . . , Nf as a quadratic function of sn (ax, k): 

<7q = ^4osn 2 (ax, fc) + Bo, g| = AjSii 2 (ax, k) + Bj. (3.5) 

Inserting ()3.5f) in (13. 4j) and equating the coefficients of equal powers of sn (ax, A;) results in the 
following relations among the solution parameters u>j, Cj, Aj and Bj and the characteristic of 
the optical lattice V$, a and k: 

A = a2k2 - mFVo , YA j = ^(±--^)-^(l- 9 ™), (3.6) 
mFffBF f-f 0BF \ m B rn F g B F/ 9BF \ 9BF J 



UJ 



3=1 

a 2 (k 2 + l) a 2 A: 2 B 

h 9BB-D0 + ^BF > Bi + — . 

2m B r-f 2mB A 



a 2 (fc 2 + l) a 2 * 2 ^- 

^ = ^2^ + ^ 0+ 2^4? (3 - 7) 

e 2 = ^(A) + B )(A + B k 2 ), e 2 = + b j)( a j + s ^ 2 )= ( 3 - 8 ) 

where j = 1, . . . , iVy. Next for convenience we introduce 
B = -A)A), Bj = -PjAj, j = l,...,N f , 
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Table 1. W = g B Fm F W B / (m B W F ) . 



1 


0o<O 


0j<O 


A > 


A,- > 


9BF ^ 


5BB S 


Vq $s a 2 k 2 /nip 


2 


/?o<0 


l<(3j< l/k 2 


A > 


Aj < 


9BF ^ 


5BB ^ 




3 


1 < 0o < 1/k 2 




A <0 


Aj > 


9BF ^ 


5BB ^ W 


Vo ^ a 2 k 2 /mF 


4 


1 < 0o < l/k 2 


l<0j< l/k 2 


A < 


Aj<0 


9BF ^ 


9bb § W 


Vq a 2 k 2 /m-p 



then 



^ = a 2 A 2 o (0 o - 1)(1 - A)fe 2 ), Cj = aSApjiPj - 1)(1 - /3^), j = l,...,N f . 



In order for our results (|3.5p to be consistent with the parametrization (|3.ip - (|3.3p we must ensure 
that both qo(x) and @o(x) are real- valued, and also qj(x) and Qj(x) are real- valued; this means 
that Cl > and gg(a:) > and also C 2 > and > (see Table H Wb = (a 2 *; 2 - m B Vb), 

Wf = (a 2 k 2 — mpVb)). An elementary analysis shows that with Z = 0, . . . , Nf one of the following 
conditions must hold 

a) Ai>0, A<0, b) Aj<0, l<A<p- 

Although our main interest is to analyze periodic solutions, note that the solutions \IH 
in (|2.1|) . (|2.2|) are not always periodic in x. Indeed, let us first calculate explicitly @q(x) and 
&j(x) by using the well known formula, see e.g. |27j : 



/ 

Jo 



du 



o P(au) - p(av) p'(av) 



1 , <r(mi — ecu) 

2x((av) + - In £ 

a a(au + cry) 



where p, £, a are standard Weierstrass functions. 

In the case a) we replace v by ivo and v by iuj, set sn 2 (iavo;k) = 0o < 0, sn 2 (iavj; k) 
f3j < and 

ei = i(2-fc 2 ), e2 = i(2fc 2 -l), e 3 = -1(1 + ^), 

and rewrite the l.h.s in terms of Jacobi elliptic functions: 

du sn 2 (iav ; k)sn 2 (au; k) 2 (in 



Iq sn 2 (iav; k) — sn 2 (oru; k) J sn 2 (cra, k) — 0o ' 

and for j = 1, . . . , Nf we have 

f x du sn 2 (iav; k)sn 2 (au; k) 2 [ X ^ U 

Jo sn 2 (iav; k) — sn 2 (au; k) 3 3 J sn 2 (au, k) — 0j 

Skipping the details we find the explicit form of 

. . _ f x du i, a(ax + iav Q ) 

&o(x) = C / . , „, 7T — = -T x + - In ■ 



o Ao(sn 2 (au;k) — 0o 2 a(ax — iavo) ' 

o / 

r = iaC(iauo) + tV-0o( 1 ~ /3o)(l - & 2 A))- 

and for j = 1, . . . , iVy we have 

A ' ~ 3 Jo ^(sn 2 (an; k) - 3 ) ~ 3 + 2 a(ax - iar;,-) ' [ ' 



Exact Solutions for Equations of Bose- Fermi Mixtures 



5 



Tj = ia((ia Vj ) + SjLJ-Ml-fyXl-tfPj). 

Pi 

These formulae provide an explicit expression for the solutions \P 6 , \£ 1 ■ with nontrivial phases; 
note that for real values of vq Qq(x), va ®j(x) are also real. Now we can find the conditions 
under which Qj(x,t) are periodic. Indeed, from (|3.9|) we can calculate the quantities To, Tj 
satisfying: 

e (x + T )-e (x) = 2tt Pq , @ j (x + T j )-@ j (x) = 2 m , j = l,...,N f . 

Then ^> b , will be periodic in x with periods To = 2mouj/a, Tj = 2rrijUj/a if there exist pairs 
of integers mo, po, and rrij, pj, such that: 

= -7T [av Q ((u) +wr /ap 1 , ^ = -Tr[avj((u) +u>T j /a}~ 1 , j = l,...,N f . 
Po Pj 

where uj (and to') are the half-periods of the Weierstrass functions. 



4 Type B nontrivial phase solutions 

For the first time solutions of this type were derived in [151 [24"1 [25] for the case of nonlin- 
ear Schrodinger equation and in |18| for the n-component CNLSE. For Bose-Fermi mixtures 
solutions of this type are possible 

• when we have two lattices Vq and Vf, 

• when m-Q = m-p. 

We seek the solutions in one of the following forms: 

q 2 , = Agsn (ax, k) + Bq, q 2 = Ajsn(ax,k) + Bj, (4.1) 
q 2 = A)Cn (ax, k) + Bo, q\ = AjCii(ax,k) + Bj, (4.2) 
ql = A dn(ax, k) + B , q\ = A,dn (ax, k) + Bj, j = l,...,N f . (4.3) 

In the first case (14. ip we have 

3a 2 k 2 3a 2 k 2 

Vb = , Vf — 

8mB 8tof 

A Q 2 ^ 2 Bj Bx _ _ B Nf S^a - — A ° 9BB 

4m F g B F Aj ' A x A Np ' 4^ j 4m B gBF A s>bf 

a 2 (k 2 + l) a 2 k 2 B 2 
wo = — - r 9bbBq + gspB 



8tob 8?TIb Aq ' 

(fc 2 + l) , a r> ^ 2 k 2 B 2 

— h 5BF-D0 - "3 72) 

8mp 8mp ^4 



6^ = ^(^o 2 " ^ " B 2 k 2 ), Q 2 = ^(B 2 - A 2 j)(A) - B)k 2 ). 

We remark that due to relations ^ = • • • = -^j- we have that all qj of the fermion fields are 
proportional to q\. 
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5 Examples of elliptic solutions 

Using the general solution equations (|3.6[) - (|3.8p we have the following special cases: (these 
solutions are possible only when we have some restrictions on g BB , 9bf, and Vq see the Table [T|) 

Example 1. Suppose that Bq = Bj = 0. Therefore we have 

qo{x) = \/~Aqsii (ax, k), qj = y/Ajsn (ax, k), (5-1) 

Ao = a2k2 - mFV °, y Aj = ^(±-^-)-^- fi-^B y (5 . 2) 

m FffBF ^ 9BF \™B m F g FB J g BF \ 9BF J 

For the frequencies ojq and u)j we have 

a 2 (l + k 2 ) a 2 (l + k 2 ) 



OJQ 



2me 2mp 
as well as Co = Cj =0. 

Example 2. Let Bq = —Aq and Bj = —Aj hold true. Then we have 



qo(x) = \/—Aocn(ax,k), qj(x) = ^/—Ajcn(ax,k). (5-3) 

The coefficients Aq and Aj have the same form as (|5,2p . The frequencies loq and ujj now look 
as follows 

a 2 (l-2k 2 ) a 2 (l-2k 2 ) 

wo = —h, " + Vq, = — ±- '- + Vq. 

2m B 2m F 

The constants Co and Cj are equal to zero again. 

Example 3. Bq = —Ao/k 2 and Bj = —Aj/k 2 . In this case we obtain 

qo(x) = — ° dn (ax, k), 1j(x) = — ^ J dn (ax, fc), 

a 2 (k 2 -2) Vq a 2 (k 2 -2) Vq 

U0= 2m B + F' Uj= 2m F W (5 ' 4) 

As before C = Qj = 0. 

Example 4. I?o = and Bj = —Aj. The result reads 

qo(x) = v^Aosn (ax, k), Qj(x) = \f— ~Ajcn (ax, k), 

a 2 (l~k 2 ) A a 2 
= — o 1- Vo + Ao9bb, Wj = ~ ■ ( 5 - 5 ) 

By analogy with the previous examples the constants Aq, Aj, Co and Cj are given by formu- 
lae (|5.2p and Co, Cj are all zero. 

Example 5. Bq = and £>j = — Aj/k 2 . Thus one gets 



qo(x) = W Aqsii (ax , k) , qj(x) = — dn(ax,k), 

k 

a 2 (k 2 -l) Vq A g BB a 2 k 2 
W ° = 2tob" T 2 k 2 ' U * = 2^ (5 ' 6) 
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Table 2. Trivial phase solutions in the generic case. We use the quantity W = 3bf™fW / b/("t-bWf)- 



1 


qo = \/A sn (ax, k) 
qj = yJ~AjSii (ax, k) 


5BF ^ 


5BB § W 


Vb 5s a 2 k 2 /m F 


2 


qo = V~ A cn (ax, k) 
qj = y 7 — Ajca. (ax, k) 


5b F ^ 


5b B § W 


Vo ^ a 2 k' 2 /m F 


3 


qo = \/— j4ocm (ax, k)/k 
qj = ^—Ajdia(ax, k)/k 


5b F ^ 


5b B § W 


Vo ^ a 2 k' 2 /m F 


4 


qo = \/A sn (ax, k) 
qj = sj — AjCQ. (ax, k) 


5BF ^ 


5b B ^ W 


V 5s a 2 k 2 /m F 


5 


qo = \A4osn (ax, k) 

qj = y/—Ajdn (ax, k)/k 


9bf ^ 


5b B ^ W 


Vo 5s a 2 k 2 /m F 


6 


qo = yj— Aoca. (ax, k) 
Qj = \/~Aj sn (ax, k) 


9bf ^ 


5B B ^ 


V ^ a 2 k 2 /m F 


7 


qo = \J— Aocn. (ax, k) 
qj = ^—Ajdn(ax,k)/k 


5b F ^ 


5b B ^ W 


Vo ^ a 2 k 2 /m F 


8 


qo = V~ Aodn (ax, k)/k 
Qj = \/Aj sn (otx, k) 


5b F ^ 


5b B ^ 


Vo ^ a 2 k 2 /m F 


9 


qo = a/— Aodn (ax, k)/k 
qj = sj — Aj-cn (ax, k) 


5b F ^ 


5B B 5s IV 


V) ^ a 2 k 2 /m F 



Example 6. Let So = ~~ A) an d -Bj = 0. Hence we have 

(/o(^c) = \/ — -^ocn (ax, fc), Qj(^) = y/Aj sn (ax, k), 

a 2 . a 2 (l-k 2 ) 
wo = 7; 5BbA), Wj = — h Vo- 

Z77T.B lm F 

Example 7. Let £?o = —Aq and £>j = —Aj/k 2 . We obtain 

qo(x) = y — Aocn (ax, k), Qj(x) = — - — dn(ax,k), 
V a 2 1 - k 2 . a 2 k 2 

W0 = 72 " 7, 1 Jn—Ao9BB, Vj = Vo - - 

k z Zm-Q k z 2m F 

Example 8. Suppose Bq = —Ao/k 2 and Bj = 0. Then 

qo(x) = ^ ° dn (ax, k), Qj(x) = \/AjSn (ax, k), 

_ a 2 k 2 gsBAp a 2 (k 2 - 1) V 

2m-B k 2 ' 3 2m F k 2 ' 

Example 9. Let B = -Ao/k 2 and Bj = —Aj. Thus 

qo(x) = ^ ° dn (ax, k), Qj(x) = \J — Ajcw (ax, k), 



a 2 k 2 k 2 -l Vo a 2 



wo = Vb - 1 r^— 5bbA), Wj = — r 



2me ' k 2 ym3JiU ' J &2 2m F 
All these cases when Vq = and j = 2 are derived for the first time in [11] . 
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5.1 Mixed trivial phase solution 

Example 10. When B = 0, B ± = 0, B 2 = -A 2 , Bj = -Aj/k 2 , j = 3, . .. ,N f the solutions 
obtain the form 

qo = \/~A~osn (ax, k), qi = \^A\sn (ax, k), 

q 2 = y / -A 2 cn(ax,k), qj = ■ s /-Ajdn(ax,k)/k, j = 3,...,Nf. 
Using equations (I3,6|) - (|3.8p we have 



9 9 N 

a z k z - Vqiuf 

Aq = , 

m F g BF 

j 



f 

E^ = ^ 2 f— 9 -^)-vo(—-¥) 

^ \m B gBF rn F g^ F J \g B F 9 BF J 



« 2 ( fc2 ~ 1) , 9BF i a i /-i \ , ffBB^O V 

^ = 2m B + + + 



« 2 (l + fc 2 ) 1 2 a 2 fc 2 

wi = — , lo 2 = a, Wj = , j = 3,...,N F . 

Imp Imp Imp 

Example 11. Let Bq = B\ = and Bj = —Aj where j = 2, . . . ,Nf. Therefore the solutions 
read 

qo(x) = \/A~Qsn(ax,k), qi(x) = \f~A\$n (ax, k), Qj( x ) = %/— Ajcn (ax, k). 
Then we obtain for frequencies the following results 

a 2 (l-k 2 ) A . a 2 (l + k 2 ) a 2 

W = |r h V +gSBA +9BFM, LOi~ 



2m B 2mp 2m F 

Example 12. Suppose Bq = —Aq, B\ = 0, B 2 = —A 2 and Bj = —Aj/k 2 where j = ?>,... , Nf. 
The solutions have the form 



qo(x) = \J — Aqcxs. (ax, k), qi(x) = yM^sn (ax, k), 

q2(x) = \J — A 2 cvl (ax, k), qj(x) = y/—Ajdn (ax, k)/k. 

The frequencies are 

V a 2 1-k 2 a\,9bf, ^ , a 2 (l-fc 2 ) 

wo = tt ~ n 1 r^(#BB-4o + 9bfA 2 ) + -prAi, u\ = Vq-\ . 

k z 2m B k z k z 2mp 

a 2 (l-2k 2 ) a 2 k 2 

u 2 = V + — , u j = V -——. 

zmp z?rip 

Example 13. Let B = -A , B x = -A x and Bj = -Aj/k 2 for j = 2, . . . , Nf. Then 



qo(x) = \/—Aocn(ax,k), qi(x) = \J — A\cu (ax, k), Qj(x) = ^—Ajdn(ax,k)/k, 
Vo a 2 1-k 2 .. a 2 (l - 2k 2 ) 

w ° = F " 2^ + HF" ^bbAo + bbfA,) , Wl = v + 2mp , 

a 2 k 2 

Example 14. Let B = -A /k 2 , B x = -A x and Bj = -Aj/k 2 for j = 2,... ,N f . Hence 

qo(x) = \/—A()dn (ax, k)/k, qi(x) = \f-A\cn. (ax, k), 1j(x) = \J — Ajdn (ax, k)/k, 

a 2 (k 2 -2) V 1-k 2 . 
w = 2mB + F + ^^(9bbA + g BP A 1 ) 1 

Vq a 2 V a 2 (k 2 -2) 



u l = To ~ 7i ) w 7 = To + 



fc 2 2mp ' J k 2 2mp 
Certainly these examples do not exhaust all possible combinations of solutions and it is easy 
to extend this list. 
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6 Vector soliton solutions 

6.1 Vector bright-bright soliton solutions 

When k — > 1, sn(ax, 1) = tanh(ax) and Bq = —Aq, Bj = —Aj we obtain that the solutions 
read 

q = V-A) T7 7, Qj = V~ A J r? V' 

cosn(ax) cosn(ax) 



where Aq < as well as Aj < 0. Using equations ()3.6j) f|3.8[) we have 

a 2 - V m F 2 

t4o = , V = Votanh (ax), 

m F g BF 

~ #BF V™B m FffBF / ffBF V #BF / 

wo = Vb - 77— a 2 , = Vb - 77— a 2 . 

IrriB Imp 

As a consequence of the restrictions on Aq and Aj one can get the following unequalities 

a 2 (a 2 - m B Vo)m F 

<?BF > 0, Vo > , f/BB < TTT 9BF, 

raF (a z — m F VQ)m B 

a 2 (a 2 - m B Vb)m F 

ffBF < 0, Vo < , 9BB > T^2 77^ 9BF- 

m F (ar — m F VQ)m B 

Vector bright soliton solution when Vo = is derived for the first time in 



6.2 Vector dark-dark soliton solutions 

When k — > 1 and Bq = Bj = are satisfied the solutions read 

qo(x) = vA) tanh(ax), = \/A* ann ( QX )- 

The natural restrictions Aq > and j4j > lead to 

^ n ^ (« 2 - m B V r )m F 2 , 

ffBF > 0, <?bb < 7^-5 777 5bf, Vo < a /m F , 

(cr — mFKoJ^B 

z n ^ (" 2 - m B Vo)m F 2 

5bf < 0, <?bb > 7^-5 777 9bf, Vo > a /m F , 

(cr — mFKoJ^B 

^F5BF ^ Sbf \ m B rn F g FB J g BF \ g BF ) 



For the frequencies uiq and u)j and the constants Co an d Cj we have 
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6.3 Vector bright-dark soliton solutions 

When k — ► 1, Bq = —Aq and Bj = 0, we have 



qo(x) = qj(x) = y r Altanh(ax), 

cosh(ax) 



a 2 



w = - 9bbA , w 3 - = V , Co = Qj = 0. 

The parameters j4o and Aj are given by (|6. 1 j) . In this case we have the following restrictions 

^ n ^ (a 2 - m B Vfo)m F T , . 2 / 

ffBF > 0, g B B > 7-5 T7~\ 9bf, V >a m F , 

(or — m-pvQ)m-Q 

^ n ^ (« 2 - m B V )m F 2 

5bf < 0, #bb < 7-5 777 9bf, V <a m F . 

(cr — mpVo)mB 

6.4 Vector dark-bright soliton solutions 

When k — > 1 and provided that Bq = and -Bj = — the result is 



I -\J ^4. j 

qo(x) = VA)tanh(ax), qj(x) = — — — -, luq = Vq + ^oflBB, 



a 



2 



cosh(ax) ' 2m F 

By analogy with the previous examples the constants Aq, Aj, Co and Cj are given by formu- 
lae (|6,ip and (|6,2p respectively. The restrictions now are 

^ n ^ (« 2 - m B Vo)m F T/ 2 / 

0BF > 0, *?bb > 7-5 777 5BF, Vb < a /m F , 

(cr — mFl/o)"iB 

z n ^ (a 2 - m B Vb)m F 2 

5BF < 0, #BB < 7— n 77T 3BF, Vb > « /"i F . 

(cr — mFl/o)"iB 
6.5 Vector dark-dark-bright soliton solutions 

Let Bq = B\ = and Bj = —Aj where j = 2, . . . , N*. Therefore the solutions read 



qo(x) = v Aq tanh(ax), qi(x) = y Ay tanh(ax), qj{x) = y— Aj sech(ax) 
Then we obtain for frequencies the following results 



u)q = Vq + qbbAq + qbfAx, iO\ 



a 2 a 2 



m F 2m F 
These examples are by no means exhaustive. 

6.6 Nontrivial phase, trigonometric limit 

In this section we consider a trap potential of the form Vtrap = Vo cos(2ax), as a model for an 
optical lattice. Our potential V is similar and differs only with additive constant. When k — > 0, 
sn (ax, 0) = sin(arr) 

(?q = Aq sin 2 (ax) + Bq, q 2 = Aj sin 2 (ax) + Bj, (6.3) 
V = V sin 2 (a2:) = -(V - V cos(2ax)), (6.4) 
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Table 3. W = g BF m F W B / (m B W F ) . 



1 


A)<o 


Pj<0 


A >0 


A > o 


5BF ^ 


SBB 5s 5BF 


^0^0 


2 


/?o<o 




A) > o 


Aj < 


5B F ^ 


9BB ^ 5BF 


Vb^O 


3 


Po > 1 


^•<0 


4) < 


Aj > 


SBF ^ 


5BB ^ 5BF 




4 


Po > 1 


/3f>l 


A) < o 


a,-<o 


9BF ^ 


9BB 5s 5BF 


V ^0 



Using equations (|3.6p - (|3.8p again we obtain the following result when (see Table [3j) 



An 



LOQ 



p2 
L 



.3. 

9bf' 



£4 



Vq L _ 5BB 
flBF V 9BF 

Nf 



a 2 + -Bo^bb + 5bf y^Bi, Uj = - a 2 + g B FB 

Im-Q z — ' zmp 



i=i 

a 2 B ( A + B ) , 6 2 = a 2 Bj (Aj + 



where 



@o(x) = arctan [ \j — — - tan(ax) j , &j{x) = arctan [ w J ^ — - tan(ax) ] . 



Bn 



B, 



>o j V V ^ 

This solution is the most important from the physical point of view [8]. 



7 Linear stability, preliminary results 

To analyze linear stability of our initial system of equations we seek solutions in the form 
ipo(x,t) = (go 0*0 + e(f>o(x,t))exp \-~^ t + i®o{x) +iK ^j , 
if)j(x,t) = {qi(x) + ecf)j(x,t))exp (—^-t + i@i(x) + itii 



and obtain the following linearized equations 
/ A Ui U 2 . . . U Nf \ 



Vi Ai 
V 2 A 2 

Vv Nf 



A Nf J 



( *o \ 



where 



A = 

5 = 

51 = 



So L - 
Lo,+ Sq 



^B(?0 V^O 



e ' -a. ' 1 



/ 

\U ,j o 







if e 2 \ 



1 



S 2 " 

~ ?f ^ ~ y ~ 35BB<7 ° ~ 9BFq ^ + 



m F qj \qj 



L, 



1 



2niF 



dl x -^j + V + g B Fq 2 ~ Wj, 



n 1 
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1 / e 2 \ 

Uij = -2g BF q qj, L j+ = \ d 2 x £ - V - g BF q^ + Uj, J = 1, ■ ■ ■ ■ 

2m F y qfi J 

The analysis of the latter matrix system is a difficult problem and only numerical simulations 
are possible. Recently a great progress was achieved for analysis of linear stability of periodic 
solutions of type (|3.1j) . (j3.2|) (see e.g. [IH |24j [25J [I8j [26] and references therein). Nevertheless 
the stability analysis is known only for solutions of type (|5.ip ~ (|5.6p and solutions with nontrivial 
phase of type (16. 3ft and (|6.4p . Linear analysis of soliton solutions is well developed, but it is out 
scope of the present paper. 

Finally we discuss three special cases: 

Case I. Let Bq = Bj = then for j = 1, . . . , Nf and qo = V^osn (ax, k), qj = yfAjsn (ax, k) 
we have the following linearized equations: 

= ~2~^x0o + M 1 ^ + 9BbAq + <7bf ^2 Aj \ sn 2 (ax,k)(/)l-uj 4> 1 , 

H\t = 2^^^ ~ + ^bb^o + 5b F ^ A, j sn 2 (ax, k)<j>$ 
+ w ^o - 2£r B FA)sn 2 (a2;, k) ^ 

3 

N>f,t = + + SbfA)) sn 2 (ax, k)4>) - ujjcfi), 

Hl t = 2^-^*0? - ( y o + 5bfA}) sn 2 (ox, k)4>J + Wj R - 2 5B F V / A)^'sn 2 (ax, k)<f$. 

Case II. Let _E?o = — Aq, Bj = —Aj then for go = V —AqCii (ax, k), qj = yJ—Ajca. (ax, k) we 
obtain the following linearized equations: 



^t = _ 2~^a!0o+ + 5bb A) + 5bf ^ A,- j sn 2 (ax,k)(f)l 



j + W I 00) 



5bbA) + 5b f y^.Aj + , 

j 

~ d ^o + ^ 3 0bb4) + 5b f + w j 



V& + 3c/bbAo + 5BF A i sn 2 (ax,k)$ + 2g BF A (l - sn 2 (ax,k))^2(f>J, 
j J j 

H4>f >t = --fy^dtxQ) + (Vo + 5BFA )sn 2 (ax, k)(j)) - (g BF A Q + u)j)(j)), 

H),t = ^^^f ~ ( y o + 5BFA )sn 2 (ax, k)<f>f + (g BF A + Uj)<j% 
- 2g BFy / A A j (1 - sn 2 (ax, k)) j = l,...,N f . 

Case III. Let Bq = —Ao/k 2 , Bj = —Aj/k 2 therefore the solutions are 



qo = a/— A)dn (ax, k)/k, qj = \J — A^dn (ax, k)/k, 
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and we obtain the following linearized equations 



= ~2m^ x< ^ )+ + 9bbA ° + 9BF A i j sn^a^Mo 
0bbA o + 5bf A i + fc2w ° j pF' 

if \ 4> R 

Hl,t = ^^ d xx0O + I 3 5BB^0 + 5BF J] A,- + /C 2 W J p-, 



[tr , q 4 , 4 1 2/ / \ iR . 2g B F^o(l ~ fc 2 sn 2 (q,fc)) \- R 
V + 3g B BA + gBF A i I sn ( ax > k )9o H p 2_> fa > 



Hit = ^ xx< pf - (V + 9 bfA ) sn 2 (az, k)<j>f + g^^+jVffi 



2 5B F y^oA,- (1 - fe 2 sn 2 (a, fc)) < 



/,., j = l,...,iV>. 

These cases are by no means exhaustive. 



8 Conclusions 

In conclusion, we have considered the mean field model for boson-fermion mixtures in optical 
lattice. Classes of quasi-periodic, periodic, elliptic solutions, and solitons have been analyzed in 
detail. These solutions can be used as initial states which can generate localized matter waves 
(solitons) through the modulational instability mechanism. This important problem is under 
consideration. 
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